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A Brief Description of the Technology 



(1 ) Any grant applications or manuscripts describing the technology: 

(4) P. K. Saha, B. R. Gomberg, and F. W. Wehrli, "A novel theory and algorithm of 
fuzzy distance transform and its applications", in Proceedings of SPIE: Medical 
Imaging, San Diego, CA, 4684, 134-145, 2002. -^J^JNe oS OoV^. , Soog^ 

(5) P. K. Saha, F. W. Wehrli, B. R. Gomberg, and M. Takahashi' Trabecular bone 
thickness from in vivo MRI using fuzzy distance transform", in Proceedings of 
International Society for Magnetic Resonance in Medicine, pp. 146, Honolulu, HI 
2002. 

(6) P. K. Saha, F. W. Wehrli, and B. R. Gomberg, "Fuzzy distance transform — 
theory, algorithms, and applications" Computer Vision and Image Understanding, 
accepted under revision. 

(2) Current curriculum vitae (CV) of inventors: Attached on separate sheets. 

(3) Closest related publications by you or others, or other prior art: 

(4) Rutvitz, D., "Data structures for operations on digital images," in Pictorial Pattern 
Recognition, pp. 105-133, 1968. 

(5) G. Levi, U. Montanari, "A gray-weighted skeleton," Information and Control, 17, 
pp. 62-91, 1970. 

(4) A concise description of the technology: 

The thickness of trabecular bone is an important architectural determinant of bone 
strength [1]. Corticosteroid exposure, for example, is well known to cause trabecular 
thinning [2]. Conversely, treatment with anabolic agents, such as parathyroid hormone, 
causes trabecular thickening. Recent advances in micro-MRI now allow visualization of 
trabecular networks in vivo. However, accurate measurement of trabecular thickness 
poses significant hurdles since in vivo resolution is usually not sufficient to fully resolve 
the structures of interest. 

Distance transform (DT) provides depth information at any point inside a structure. The 
method described in this invention is based on sampling DT values along the medial axis 
thus yielding a regional thickness measure. However, whereas DT has previously been 
applied to binary images (see [3, 4] for early intensity based DT) such representations 
fail to handle object material heterogeneity and partial volume effects. Partial volume 
blurring is particularly severe at in vivo resolution of MRI and CT images of trabecular 
bones since the voxel size is typically comparable to or larger than trabecular thickness. It 
is well known that under these conditions the images cannot be binarized and the DT of 
hard objects fails (see, for example, Wehrli FW, Hwang SN, Song HK, Gomberg BR 2000 
Visualization and analysis of trabecular bone architecture in the limited spatial resolution 
regime of in vivo Micro-MRI. In: Maiumdar S, Bay BK (eds.) Noninvasive assessment of 
trabecular bone architecture and the competence of bone. Kluwer, Big Sur, California, pp 
153-164). The present invention overcomes the limitations intrinsic to binary distance 
transforms with a method capable of generating distance transforms of fuzzy objects 
referred to as fuzzy distance transform (FDT). The method can be understood by 
considering the paths between two points in a fuzzy object. A path n from a point p to 
another point q (not necessarily distinct) is a sequence of (p = p,,^,...,p ( = q) of points 
such that every two successive points in the sequence are adjacent (spatially nearby) to 
each other. The length of n is the minimum amount of material needed to traverse to 
complete a walk along the path. It is calculated as follows: 



length of n = Xl(Mp,) + /Up /+1 ))x|| p,-p M \\ , 



where \i a is the material density (here, bone volume fraction, short, bvf) function and [|«|| 
gives the Euclidean distance. There are infinitely many paths between two points. The 
fuzzy distance between two points is the length of the shortest path between them, i.e. the 
one requiring the minimum amount of material to be traversed from one point to the other. 
It has been shown in [5] that fuzzy distance satisfies the desired metric property. 

The support of an object is the set of points with nonzero material density (here, bvf). FDT 
is an operation that assigns at every point its fuzzy distance from the background - the 
inverse of the support. The method has been reduced to practice and validated in 
trabecular bone images of varying resolutions, orientations, and signal-to-noise ratios. 
Specifically, a dynamic programming-based algorithm for efficient computation of FDT 
has been conceived and implemented. Let Q. 0 (p) denote the FDT value at p in a fuzzy 
object (here, bvf mapped image of TB) 0 and let Sk{@{0)) denote the skeleton of the 
support of O . Sk(®{0)) can be computed using currently existing methods. The 
average thickness of O , denoted by x{0) is then computed as 
£ 2Q*(P) 



T(<?) 



\sk(e(&))\ 

where, | Sk(Q(6>)) | yields the number of points in the skeleton. Although, the target 
application is trabecular bone images, the method is applicable to images of wide range 
of elongated structures possessing heterogeneous material distributions and/or fuzzy 
boundaries generated either by true object properties or by partial volume blurring due to 
the limited resolution of imaging device. 
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